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Abstract: This paper proposes a multiple-objective fuzzy 
optimization model for Independent Transmission Company 
(ITC) to plan its transmission enhancements. There are 
several objectives to be considered, such as investment cost 
objective, power transferring efficiency objective, revenue 
objective, network congestion mitigating objective and 
environmental concern objective. The uncertain capital 
investment are also taken into account in the proposed 
model. By applying the fuzzy set theory, the above multi- 
objective optimization model has been converted into regular 
single-objective optimization problem to be solved. The 
numerical examples are given to demonstrate the 
effectiveness of the proposed model. 
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I. INTRODUCTION 


Transmission plays a key role in making the competitive 
electricity market work. It provides the means to broaden 
and strengthen competitive generation markets. A more 
robust transmission system would not only bring in 
competitive bidders from far away, but would also 
debalkanize market pockets in which dominant generators 
can exercise market power due to transmission restraints. 
However, transmission has not kept up with the growth in 
market, and some of the favored organizational fixes seem to 
do no more than ration existing assets " I, Creating the 
necessary transmission organizations and providing them 
with incentives to meet the demands placed on them is vital 
for the robust competition. The Independent transmission 
company (ITC) has offered a best structural form for 
transmission. 


ITC is a commercial, profit-making transmission entity. It is 
independent from generation ownership and control. ITC 
owns and operates transmission assets, and it has rights to 
plan, finance and execute transmission expansions. ITC can 
benefits from increased efficiency or growth of market, and 
suffers commercial penalties for inadequate reliability or 
service. ITC has the ownership and operation together, and 
this efficient business structure encourages optimal decision- 
making. With the proper pricing framework and regulations, 
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ITC can expand the network to carry more volume, 
and reduce congestion, and then offers an efficient, 
economical and reliable transmission for electricity market. 


Although many approaches °4*) have been proposed to plan 
the transmission networks in traditional vertically regulated 
power industry structure, it is still necessary to study how to 
plan the transmission expansion of ITC in the market 
environment. 


This paper proposes a multiple-objective fuzzy optimization 
model for ITC to make its transmission enhance decisions. 
The maximization of network congestion mitigating 
objective, the maximization of power transferring efficiency 
objective, and the minimization of environmental impact 
objective are proposed to constraint ITC to offer good 
services to market customers and less environmental 
impacts. The minimization of capital investment, the 
maximization of revenue objective are considered for its own 


commercial success of ITC. The risk of capital raising of 


ITC is taken as constraint to be considered. By employing 
the fuzzy set theory, this multi-objective fuzzy 
optimization problem has been converted into regular single- 
objective optimization problem to be solved. The numerical 
examples are given to demonstrate the effectiveness of the 
proposed model. 


Il. INVESTMENT DECISION-MAKING 


The goal of transmission enhancement planning is to 
determine to add or enhance what type of lines in which 
routes of power systems. 


It is assumed that all possible enhanced or new-added routes, 
and all possible candidate line types in each route are given. 
The line types can be varied with line rated voltage, 
conductor size and type, tower configuration, and different 
parallel circuits. 


- Based on DC power flow models, the following multi- 


objective fuzzy optimization model is proposed to make 
expansion decisions. 


1.Objective functions 
(1). The first objective is to minimize the investment costs: 


Min Z, =C(X) (1) 
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where, C(X) = > YC, X y Tepresents the investment cost 
ieM jel, 

of expansion plan. X stands for the set of investment 

decisions. C,, is the investment cost of candidate line type j 


in route i. X, is 0-1 integer decision variables. M is the 
set of routes. Z; is the set of candidate line types in route / 


(2). The second objective is to maximize the power 
transferring efficiency. The minimization of transmission 
loss is used to represent this objective: 


Min Z, = D(X) (2) 
where, D(X) is the expected power loss of expansion plan. It 
is approximately expressed as the linear function of 
investment decisions X as follows: 


fe)8) = = aD 
DXY)=D,+ > = PiXj-Pi [+ > —— BX, —B, | 3) 
By pis [5 my Lay [s te 


where, D, is the expected power loss of existing network, 
aD oD 
~—== and —— are 
OP, i OB; 


respect to the active capacity P; and susceptance 8, of 


the sensitivity factors of power loss with 


route i in fhe existing network, Py and B, are, the 
active capacity and susceptance of candidate line type j in 
route . 


(3). The third objective is to maximize mitigating the 
network congestion. It is expressed as the minimization of 
system expected demands not served(EDNS): 


Min Z, = E(X) (4) 
where, E(X) is the system expected demand not served, and 


it is approximately expressed as the linear function of 
investment decisions X : 


OE 
— > By Xz - B, 35 
aB, [> FO i I ) 


aE { o> = 
E(X) = Ey + Y= | Py X; ~ Pi |+ 
ieMOP' \ jal, ieM 


where, E, is the EDNS index of existing network, = and 
i 


&. are the sensitivity factors of EDNS with respect to the 
i 
capacity and susceptance of route 7 in the existing network. 


(4). The fourth objective is to minimize the environmental 


impacts. Since the right-of-way occupation is the main 
public concern of transmission, the minimization of line 
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right-of-way occupied area is used to express this objective. 


Min Z, = F(X) (6) 


where, F(X) = > AX y is the new-added right-of-way 
ieM jel, 

occupied area of expansion plan, F,, is new added right-of- 

way occupied areas of candidate line type j in route i. 


(5). The fifth objective is to maximize the revenue from 
transmission. Because the transmission revenue is heavily 
relied on the volume of transmission, the objective is 
expressed as the maximization of load carrying capacity of 
network: 

Max Z,; = G(X) (7) 
where, G(X) represents the load carrying capacity of 
expansion plan, it is approximately expressed as the linear 
function of investment decisions: 


GX)=G, “32 5h P) 32 Sy -a| (8) 


where, Go is the load carrying capacity of existing network, 


a and oe are the sensitivity factors of load carrying 


OP; OB; 
capacity with respect to the capacity and susceptance of 
route j in the existing network. 


2.The constraints 


(1). Capital investment limits: The expansion decisions is 
determined by how much money ITC can raised. The 
financial risk is always existing. So, the following fuzzy 
constraint is used to express the uncertain of capital 
investment: 


c(x)sc (9) 


where, C is the estimated maximum investment limit. “<” 


is a fuzzy unequal symbol, and it means “roughly less and 
equal than”, 


(2). The maximum acceptable power loss. D, limits: 
D(X)<D (10) 


(3). The maximum acceptable value of system expected 
demand not served, E , limits: 


E(X)<E (11) 


(4). The minimum. acceptable value of new-added right-of- 
way occupied area, F’, limits: 


F(X)sF (12) 


(5). The minimum acceptable load carrying capacity,G, 
limits: 
G(X) 2G (13) 


Besides above constrains, the investment decision variables 
are subject to the following constraints: 


VX, =1  ieM (14) 
jel; 
X,={0R jel), ieM (15) 
Ill. OPERATION SIMULATION 
In above model, D(X),E(X)andG(X)are directly 


associated with system operation decisions. So, it is required 
to simulate the network configurations under various system 
operation states. 


1. System state selection 


System operation states are determined by the load levels, 
generation unit commitment and line states. The total 
number of system states for a practical system is so large that 
it is impossible to simulate all the system states. So, only the 
selected typical system states are used in this paper. 


For load demands, several typical load levels are used to 
represent load variations. 


For generation units, it is associated with load levels to be 
considered. For each typical load level, one generation unit 
commitment is given correspondingly. The generation 
outages are neglected. 


For lines, both normal and outage states are considered, but 
the multiple outages is only considered up to 2 order. 


2. System operation behavior simulation 
For any system state +, the following double-objective 


optimization model is used to simulate the system operation 
behaviors: 


Min W, = >) R; (16.a) 
ieN 
Min W, = }\7,P? (16.b) 
ieM 
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st. Pg +R-BO=P, (16.c) 
P|sP (16.4) 

Pg SP <Pg (16.e) 
O<R<P, (16.f) 


where, Nis the set of system buses, 7, and P; are the 
resistance and active power of route i, Pp and R are the 
bus load demand vector and curtailment vector, Pg, Py and 
P, are the bus generation vector and upper and lower 
capacity limit vectors, B and @ are the bus susceptance 


matrix and voltage angle vector, P is the route capacity 
limit vector. 


The model is solved in two steps: 


Step 1: Only objective W, is to be optimized under the 
constraints of Eqs. (16.c)~(16.f), then the minimal value of 


system demand not served, 7," , is obtained. 


Step 2: Only objective W, is to be optimized under the 
constraints of Eqs. (16.c)}~(16.f) and W,<W,, then the 
minimal power loss, W, , is obtained. 


Ww, and W, are the system demand not served £,, and 
power loss D, in-system state ¢. 


The sensitivity factors of E, and D, with respect to route i 
capacity and susceptance can be calculated as follows: 


GE, [OPi =a5 (17) 
aE, /@B; = (ap, -&p,, (Gn — 9:2) (18) 
aD, /AP: = Bs (19) 
aD, /@B; =(Bp,, - Bo, Xn — 92) (20) 


where, as and Bs, are the dual variables of constraint i in 
Eqs. (16.d) with respect to objective W, and W, 
respectively. @p,(@p,,) and Bp, (Ap,,) are the dual 
variables of constraint 7, (i.) in Eqs. (16.c) with respect to 
objective W, and W, respectively. i, and i, are the 
terminal buses of route 7. 


3. System load carrying capacity simulation 


For any system state ¢, the following optimization model is 
used to simulate the system load carrying capacity: 


Max W; = >” Ps (21.2) 
; ieN 
s.t. Po -B@=P, (21.b) 
[P|<P (21.¢) 
Pg $Pg Pg (21.d) 
Py SOP ieN (21.e) 
where, P;,is load demand that bus i can supply, 


6;,=P,,/ X Po, »Pp is forecasted load demand of bus 7. 


ie N 
It is assumed that the proportion of every bus load to total 
system load is unchanged when loads is increasing. Pe is 


the maximum load carrying capacity under normal system 
state. P_ is obtained by solving the model (21), but Eqs. 


(21.e) is changed as an equal constraint. 


The optimal objective of the above model, W;, is G,, the 
system load carrying capacity in state ¢. And, the sensitivity 


factors of G, with respect to route i capacity and 
susceptance can be computed as follows: 
8G, /OP: =r, (22) 


aG, / OB; = (Yd, -Y Dry M6 - 42) (23) 
where, 75 is the dual variable of constraint 7 in Eqs. (21.c), 


Yp, (¥p,, ) is the dual variables of constraint i, (i, ) in Eqs. 
(21.b). 


4. The sensitivity factor calculations 


By applying above two models to all selected system states, 
we can obtain the required constants and sensitivity factors 
that used in the functions of D(X), E(X) and G(X) as 
follows: 


Dy = = B,D, (24) 
tef 
oD oD, 
OP; de OP; on 
oD. OD, 
—= —. 26 
aB, x a) ue 
Ey = >) PAE; (27). 
ief 
OE yee (28) 
OP; eT OP, é 
OE OE, 
pore’ Py, (29) 
OB, 2 OB; 
Gp = >Y 2G, (30) 
tel 
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8G aG, 


OP; or OP; oy) 
aG aG, 

— = — 32 
OB, dP «OB, seu 


where,7 is the set of selected system states, p,is the 
computed probability of state ¢. For system state f, its 
computed probability is: 
Pi =P, (3) 
yr 
teT 
p= [[a [][a-9) 
ieMc 1eM,-Mo 
where, q,is the outage probability of line 7.Mé, is the set 
of outage lines in state t. M, is the set of total lines. 9, is 
the occurrence probability of the load level of state ¢. 


(34) 


IV. SOLUTION ALGORITHM 


The fuzzy set theory is employed to solve above-mentioned 
multi-objective fuzzy optimization problem. 


First, we define fuzzy set F to express the satisfying degree 
of constraint (9). Its membership function is defined as 
follows: 


: Pe C(X)2C+c, 
pg(X)={¥S——®) Gcce<é+c, G5) 
‘4 crys 


where, the uncertain investment limit is varied within the 
range of |C,C+C,|,C, is the tolerance of investment limit. 
4g (X)is the possibility of future capital investment to be 
raised. 


Then, the fuzzy sets C,D,E,F and G are introduced to 
express the realization degrees of objective (1), (2), (4), (6) 
and (7) respectively. The membership functions of these 


fuzzy sets are defined as follows: 
a 0 C(X)2C+C, 
ue(X) = eee Cc’ <C(X)<€+¢, (36) 
ee c(x)sc" 
5 0 . D(X)2D 
#5 (X) = Ul : ) D< D(%) < D G7) 
1. D(X)s 0" 
; 2 “ E(X)2E 
p(X) =} © FeEX<E 67 
1 E(X)s E* 


0 F(X)2F 

yg) = {= = F< F(X)<F (38) 
1 F(X) Ss F" 
0 G(X) SG 

He (X) = 2” 2 G<G(x)<G’ (39) 
= G(x)sG 


1 


where, C’,D°,E',F’and G" are the optimal objective 
values of objective (1), (2), (4), (6) and (7) under the 


constraints of Eqs. (10)~(15) and C(X)<C+C,. 


Since the realization of objective (1), (2), (4), (6), (7) and 
satisfying of constraint (9) are contradictive éach other, it is 
impossible to completely satisfy all these objective and 
constraint simultaneously. So, it is the only feasible way to 
make a compromise that has the maximum overall satisfying 
degree of all objectives and constraints. 


Here, fuzzy set 4 is used to express the overall satisfying 
degree of the realization of objective (1), (2), (4), (6); (7) and 
satisfying of the constraint (9), and its membership function 
is defined as follows: 


~ mE Ew DYN 


where, sea oz is forradused, to scien the relative 
importance of abiectve or constraint , and @, 21. 


Hence, the proposed model. described in Eqs. (1}~(15) is 
equivalently converted to maximize y,(X) subject to 


mentioned routes, all other routes can add Type B 220kV 
lines, the investment is 25000$/km, the line capacity is 200 
MW. The total parallel circuits cannot excess 2 circuits per 
route in all routes, except route 4-24 and 4-5 can have 3 
parallel circuits. The line outage rate of 500kV and 220kV 
lines are p =0.0025L +0.126 (1/yr) and 
p=0.0034L +0.290(1/yr) respectively, L is the length of 
line in km, the repair duration of line outage is 16h. The 
widths of right-of-way occupied area are 48.27m for single 
220kV line, 67.58m for single 500kV line, 68.27m for 
double 220kV line, 97.58m for double 500kV line, 88.27m 
for triple 220kV line, 117.58m for triple SOOkV line. The 
following constants are used in the optimization process: 


€=$35 million, C,=$35 million, D=55.0MW, 
E=3.57MW, F =70km’, G =3410MW. 


Fig.1. The 28-bus power system 
There are 2 typical load levels to be considered, and the 
yearly equivalent duration of each level is 3500h. Table 1 
gives the bus load demands of different load levels, and 


Table 2 lists the generation lower and upper limits for each 


Eqs. (10)~(15), that is: load levels. 


Max y5(X)=4 Table 1. Busload demands ( 
Spee ef Bus [Typical fond states] 1 Teale 
ne, Ug (X)2 og ~ Kee, DERG aH ie Typical load states PB Yypical load states 
O<sA<l a at 90 a a 
Eqs. (10)~(15) C2 ao BB 300 320 
(41) I 
Then, we can obtain the optimal solution as X,( je ,, Se et SSE a Ce 
* * 6 200 ¢ 
lal) and 2, X3(Jely, 1004) te fry opis SPE tt 
satiety Grates of Hila expansion plan: 10-299 99} tr 
IO ae DC | Se See 


_Table 2. Bus generation limits (M 
eneration Limits 
(Load state 1) 


V. NUMERICAL EXAMPLES 


te 


A 28-bus system is used to demonstrate the validity of the 
proposed method. 


Fig.1 depicts its network configurations. In Fig.1, the solid 
lines stand for existing routes, the dash lines are new-added 
routes, the data above the line is the route length(km). In 


route 12-24, 12-27 and 24-28, 500kV lines can be added or able 3 end’ Table-4 pive the planning results of the system: 
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enhanced, the investment cost is 65000S/km, the line 
-capacity is 1SOOMW. In route 7-27,12-23,21-22,21-23 and 
23-24, Type A 220kV lines can be added, the investment is 
36000$/km, and the line capacity is 300 MW. Except above- 
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Plan 1~5 are the optimization results of single objective (1), 
(2), (4), (6) and (7), and plan 6 is the optimization results of 
all objectives together. Table 3 lists some characteristic 


attributes of expansion plans. Table 4 lists the new-added of one particular objective is maximized. For plan 6, the 

lines of expansion plans. satisfying degree of single objective is not the best, but the 
overall satisfying degree is the best. Plan 6 is the 

From Table 3, we see that the overall satisfying degree of | recommenced optimal expansion plan. 

plan 1~5 is not good enough, although the satisfying degree 


Table 3. Comparisons of different expansion plans 


a a a 
She en | 
32.14 U.S | RSS| 4 49.68 

Power Loss (MW) [332735453552 SiR 


DNS (MW) 2 LS 
RightoF-way oscupled area 386 | 06ST 
[Load carrying capacities (MW) | __ 3487 __|__ 4364 |~_ #366 [3486 [36 3939] 


Objective 
realization degree 


[__ObjectiveC [10 | _ 0.0046 | 0.0284 [1.0 01a 0.5374 


Objective E J 0.0290 fT 00.0063 0.9508 0.6616 
|____ Objective F—__ | 0.9369} _0:0286 | 0.050510 | __ 0.1610 | 0.6737 
uture investment possible realization degree | ss 1.0 | 0.0050 ~~ [_—_—0~.030 : 


Enhanced Routes 
220 Lines Ok’ 


jingle Circuit Dubie Circuit Single Circuit” [Double Circuit | 
0-27, 12-13, 12-43, 16-27 4-19, 19-25,21-23 P12-27,24-28 
q,5-15,5-00,6-16,9-18, 16-27 3-05,4-5,4-19,4-24,6-15,10-27,12-13,12-23, 15-22 
18-23,19-25,21-23 


3 
-2 


9-18,10-27, 12 
12-23, 16-27,21-22 
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